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ABSTRACT 

We have discussed non-linear stability in photogravitational non-planar re- 
stricted three body problem with oblate smaller primary. By photogravitational 
we mean that both primaries are radiating. We normalised the Hamiltonian using 
Lie transform as in Coppola and Rand (1989). We transformed the system into 
Birkhoff 's normal form. Lie transforms reduce the system to an equivalent sim- 
pler system which is immediately solvable. Applying Arnolds theorem, we have 
found non-linear stability criteria. We conclude that Lg is stable. We plotted 
graphs for {uJi,D2). They are rectangular hyperbola. 

Subject headings: Non-linear stability:Photogravitational:Non-planar:Oblate pri- 
mary: RTBP 

1. Introduction 

G. Hori(1966,1967) applied a theorem by Lie in canonical transformation to the the- 
ory of general perturbations. Theorem is applicable to such cases where the undisturbed 
portion of Hamiltonian depends on angular variable as well as momentum variables. A. 
Deprit(1969) introduced the concept of Lie series to the cases where the generating function 
itself depends explicitly on the small parameter. Lie transforms define naturally a class of 
canonical mappings in the form of power series in the small parameter. They reviewed how 
a Lie series defines a canonical mapping as a formal power series of a small parameter e, 
provided the generating function itself does not depend upon e. This restriction is overcome 
by introducing Lie transform. They showed that how they naturally define the canonical 
transformations contemplated by Von Zeipel's method. Canonical mappings defined by Lie 
transforms as formal power series of a small parameter constitute the natural ingredient of 
Transformation Theory applied to Hamiltonian systems. 
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Orbital stability of quasi-periodic motions in multidimensional Hamiltonian systems 
was studied by Sokolskii(1978). With some applications to the Birkhoff's normal form along 
with its generalized form by K. R. Meyer, the restricted problem of three bodies near L4, the 
Birkhoff's normalization procedure, and the singular perturbation, of Hamiltonian systems 
have been discussed by Liu(1985). K.R. Meyer and D.S.Schmidt(1986) established the full 
stability of Lagrange equilibrium point in the planar restricted three body problem even 
in the case when /i — fj,c . Hamiltonian is normalised up to order six and then KAM 
theory is applied. This establishes the stability of the equilibrium in degenerate case. A. P. 
Markeev(1966) and K.T. Alfriend(1970,1971) have shown that L4 is unstable when the mass 
ratio is equal to 112 or /X3. The Lie transform method is an efficient perturbation scheme 
which explicitly generates the functional form of the reduced Hamiltonian under an implicitly 
defined canonical periodic near identity - transformation. 

Coppola & Rand(1989) applied a method of Lie Transforms, a perturbation method 
for differential equations to a general class of Hamiltonian systems using computer algebra. 
They developed explicit formulas for transforming the system into Birkhoff normal form. 
They formed explicit nonlinear stability criteria solely in terms of H for systems where the 
linear stability is inconclusive. They applied these results to the non-linear stability of L4 
in the Circular RTBP. At L4 , Arnold's theorem(1961) must be used since a Lyapunov 
function cannot be found. They conffimed the previous computations of Deprit and Deprit 
Bartholome (1967), Meyer and Schmidt (1986). Algorithms of linear and nonlinear normal- 
ization of a Hamiltonian system near an equilibrium point were described by Maciejewski, 
and Gozdziewski(1991). 

A. Jorba(1997) described the effective computation of normal forms, centre manifolds 
and first integrals in Hamiltonian Mechanics. These kind of conclusions are very useful. They 
allow, for example, to give explicit estimates on the diffusion time or to compute invariant 
tori. Their approach presented here is based on using algebraic manipulation for the formal 
series but taking numerical co-efficients for them. 

J. Palacian & P. Yanguas (2000) described the reduction of perturbed Hamiltonian sys- 
tems. They used a technique based on Lie-transformations. By extending an integral of the 
unperturbed part to the whole transformed system up to a certain order of approximation, 
the number of degrees of freedom of such a system is reduced, under certain conditions. The 
idea of reducing a perturbed system is valid not only for Hamiltonians, but also for any 
system of differential equations. Recently there has been a resurgence of this subject. The 
problem of building formal integrals for Hamiltonian systems has received a wide treatment 
in the last forty five years. The results have been applied in fields such as Molecular Physics 
or Astrodynamics. For example, in galactic models of three degrees of freedom, the search for 
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third integral becomes very important to analyze the onset of Chaos. Their approach con- 
sists in generalising the concept of normal forms by selecting a function G{x) and proposing 
thereafter a symplectic change of variables. 

The nonlinear stability of triangular equilibrium points was studied by Kushvah et al 
(2007) in the generalised photogravitational restricted three body problem with Poynting- 
Robertson drag. They have performed first and second order normalization of the Hamil- 
tonian of the problem and applied KAM theorem to examine the condition of non-linear 
stability. After computation they have found three critical mass ratios and concluded that 
triangular points are stable in the nonlinear sense except three critical mass ratios at which 
KAM theorem fails. 

Hence, we thought to examine the Non-linear stability of Lq, equilibrium point of Non- 
planar photogravitational restricted three body problem with oblate smaller primary. Wc 
examined the linear stability of above problem in Shankaran(2011) where qi is radiation 
pressure of bigger primary and q2 that of smaller. We have found that Lq is unstable due to 
positive real part in complex roots. Now, we proceed to normalize the Hamiltonian of the 
problem as in Coppola and Rand (1989) using Lie transform. We transformed the system into 
Birkhoff 's normal form. Using Arnold's theorem, we have found non-linear stability criteria 
in terms of H. Lie transforms reduce the system to an equivalent simpler system which is 
immediately solvable. Finally we find that D2 7^ 0. Hence, according to Arnold theorem, we 
conclude that Lq is stable. We plotted graphs for (wi, D2). They are rectangular hyperbola. 



2. Computation of D2 

We suppose that q is the radiation co-efficient of bigger primary and Q that of smaller 
primary in photogravitational non-planar restricted three body problem. We want to apply 
Arnold's theorem(1961). We applied method of Lie transform using computer algebra as 
given in research article of Coppola and Rand(1989). Wc suppose L„ = {,Wn} and S'„ = 
n Sm=o -^n-m'S'mj (where 71 = 1, 2, 3 ... ) be the operators. Then we transformed the near- 
identity transformation {x^, Zm) to (X„, Zm) variables. The Kamiltonian is given by 



Hn + -{Hn^Wn} + -"^[Ln-mKrn + rnSn-mHrn] 71 = 2,3,4... (l) 

We choose the generating function Wn in such a way so as to best simplify Kamiltonian and 
each term to be canceled will be of the form AXlYlX^Y^ where A is a constant. Wc choose 
Wn to be a sum of terms, one for each term to be canceled, of the form Wn = AXlYlXl^Y^ 
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where B 



with n — j-\-l-\-r-\-s — 2. Here uii and ui^ are the basic frequencies 



uji{l—j)+u>3{s-r) 

which are rationally independent and we suppose that the frequency in vertical direction 
is constant i.e. cua = 1. If j = / and s = r then B will be infinite hence the terms 
{XiYiy {X2Y2Y . The D2 is given as (please see Coppola and Rand(1989)) 



D2 = -{K22mujl + KiniuJiU3 + K0022ojI) 
We performed computation using Mathematica and found the following results: 



(2) 



K0022 = H0022 + i 



1 1 
— HOlll X HlOll H (H0012 x H0021 



+ H0003 X H0030) 
1 



+ 



H0120 X H1002 + 



^— — H0102 X H1020 1 ; (3) 



Kllll =Hllll + i 



(HlOll X H1200 



-H0120 X H1002+ 



+ HOlll X H2100) + , . 

[uji - 2UJ3) 

— (HllOl X H0021 + HlllO * H0012) + 



H0210 X H2001 + 



1 



H0102 X H1020 + 



1 



-H0201 X H2010 1 
(2a;i + U3) J 



(4) 



K2200 = H2200 + i ( — HllOl x HlllO + — (H1200 x H2100 

\UJ3 Ui 



(5) 



+H0300 X H3000) + 



^ -H0210 X H2001 + ^ -H0201 x H2010 ) 

(2a;i-a;3) (2a;i+a;3) / 



K2200 
Kllll 



-863 + 



12aia3 , ai , ai , 60204 



(6) 



+ 



+ 



+ 



cui — 2a;3 2uji — uj^ oj^ 2uji + uj^ uji + 2003 



(7) 



K0022 = - r-1265 + ^ 
8 V i^a 



+ «3 



+ 



2uJi 



(8) 
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^2 = - (K2200a;3^ + KllllwiWg + K0022a;i2) (9) 



where 



c=(3*A)V2_£((lzgM!; 



3/2 



((-1 + //)V(-l + /^)V) + (24 (v^g//^ - v^gV' + v^QV(-l + /i)2v^- 



VSgQ^ FTTT^V ^ - 3V3Q v^FTT^v^ + 

(945 (g^u^ + Qv/(-l + /x)2^ - 6gV(-l + /x)2/xv^ - 20g ^(-1 + /i) V 
15Q + /i) Vv^ - 6Q + /i) V + 

g + /i) Vv^ + isgvi-i + / 

(S ((-1 + /.)V(-1 + /^)V)) + 



02 = 

135 (V3g 



+ 



2((-l+A')V(-l+^')'j 



40g ^(-l + ^")2 /iv/;;i2 ^ 39gg^( 1 ^ ^j fij ^^ 

40Qv/(-l + /x) 2^3^ ^ 34gQ^(-l + 10O^(-l + /^)2 

/i^ - 6gQ v^(-l + /i)2/i^ - qQ ^(-1 + /i)2/i5 ^+ 

60g^(-l+;,)2(;,2)3/2 _ 565Q^(_l+;,)2 (^2)3/2^ ^2^ 
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(g(-i + A.)V(-i + /^)V) + o[Afi^- 

n - ( 39(1-m) _ 39(1-m)m _ 3Qv^\ 

/ 45g(l-M) 452 |_ 45g(l-M)M j 45g/x , 

V 2((1-M)2)^/^ 4(1-M)({1-A')=')'/' 2((l-M)2)^/2 ^ 4(l-^)((l-M)2)5/2^ 

45QV^ ^ >1 + 36^/3(l-g)q(l-/x) _ 36^/3(l-q)g(l-M) 



36^3-^ 36^g.^ 36^.^ 36x/3q^\ .3/3 
f 945(j _|_ 945ij _ 945gM 



( 3^% _ 3V3gv/(l-;<)-^M 3v^0v^\ /T I 
\^2((1-/.)2)V2 2(l-^.)« + 2i? yV/l+ 



16(1-m)3((1-m)2) 
04 



75v/3gA3/2 j_ l 3 ^ I ^'ItT 90(l-g) 



^(-1+A')V(-1+/^)' ^ \((1-M)^)''/' Q((l-M)^)'^'/V J ^ 

135gv/(l-M)2 _ 243g2 ^(1-^)2 ^ 27g2 ^(1-^)2^ 135^ 

Q(l-M)6 2Q(l-^)6 2Q(1-M)6 + 

135q^/J^ 135^ , 243g^ , 27q^ ^ ^2 , /I15/2. 

j;^ ^ 2;[i5 ^ 2^5 I ^ +^L^J ' 

Qv/(-l + - 3QV(-1 + (/i^)'/') / ((-l + /x)V(-l + /^)V) - 

(15 (-3g/i^ - 2Q^(-l + ^)2^ + 10Qy/{-l + fiYfi^ + 20gV(-l + /^)' 



/x^^ - iooV(-i + /")V\/^ + 2QV(-l + /i)Vv^- 

2ogv(-i + /^)M/^')'^') ^) / (2 ((-i + /^)V(-i + /^)V)) + 

(30 (^^//S - V3gV' - V3gV(-i + //)2v^ + y3?gv(-i + /i)2v^+ 



^g v/(-i + /i)V + Vsgg v/(-i + /x) V v^- 

^1=^/') / (g(-i + //) v(-i + v) - 

(945 (g/x^ - Q^{-1 + + 7gv/(-l + /i) + 35gV(-l + ^if^i'^- 

35Q V(-i + /i)V + 2ig^(-i + /i)V Va^- 
7gv(-i + /^)W^ + gV(-i + /^)>W - 2igv(-i + /^)M/^'f'^ 
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A^) I (4 ((-1 + /.)V(-1 + /^)V) ) + 0{Afl\ 



3g I 3gM 



((i_^)2)6/2 -r ((i_^)2)6/2 ^ 4^ 4^7? r I -r 

945g 45g 45g^ (1-^)^ _ 945g/i , 

8((l-//)2)''/2 8(l-At)2{(l-At)2)"/2 4{1-m)« 8((l-/i)2)'^/2-r 



45gjU 



4(1-m)8 2//7 I ^ 



8(l-,,)2((l_^)2)S/2 

(^90 (^g/i^ - \/3gV^ - Vsg/FrTTO" + \/3gQV(-l + /x)2v^+ 



(4725 (^qfi^ - Q^i-i + fiy^ + 7gv/(-i + /i)>\/^ + 35gV(-i + /^)V\/^ 

35g v/(-i + /i)V + 2igv/(-i + /u)V v^- 

7gv(-i + /^)Vv^+ gyFT+TlF)" W - 2ig^/Fi+7^ (/^')'^') 

/ (4(-l + /.)V(-l + /^)V) + 

3? 3gn : SQ^/J^ \ f 453 45g^(l-M)^ , 

~r Q.,5 "T NO.., NON5/2 Afl ..\8 ~r 



8((l-/i)2)5/2 8((1-m)2)^/2 SA'^ / V 16(1-m)2{{1-/")2)^/^ 4(1-m)8 

45g/t 



16{l-/.)2((l-/,)2)5/2 



45g^(l-M)"/x 75Qy5!^ J 1 
^ 4(1-At)« 8/i'^ J 



45V3(l-g)g 45V3(l-g)g , i^VSy/Ji^ _ 45V3gy^ 

4Q((l-/i)2)6/2 4Q((l-/i)2)6/2,i 4m'' 4/x7 



4573 45v;3gv^ \ ^3/2 



4//6 + 4/i6 

945g I 315g^ 945g/t 

64(1-m)4((1-m)2)5/2 2(1-m)10 64(l-/i)4((l-/,)2)5/2 

315gV(l-M)2M 11025Qy^\ .2 , nrjl5/2. 
2(1-^)10 64/x9 I ^ + ^^[^J , 

Prom above results D2 upto first order in A is 

„ _ 9gy/il-fiy^uj'i _ 9gV(l-M)^M^f , 9Qy/^i4 

-^2 16(l-^t)6^_ 16(l-iu)6 16^5 

9gV(l-/^)V(-l+/^)^^3 9Q2a;3 _ 9Q2t^3 9Q2^3 

4(1-M)« ( -1+^IP~ ^ 4(-1+ m)2/.6 2(-1+^2^5 4(_i+/,)2^4- 

2(1-/.)«(-1+m)^ + 4(1-m)^(-1+/x)4 
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(1 -/^)6 (-1 +^)^3 2(l-At )6(-l+M)V2 + (l_^)6(_l+^)2^ 



2(-l+At)-''' 2(-l+/i)V''' ^ 2(-l+At)3M'' 2(-l+/i)'V3 

2(-l+;t):V^ iPTTTI)^ 4(-l+/i)4^6tJi + (-l+/j)4^ 5a;7 



{-I+m)^/^^^! 2(-l+At)V2a;i 16(1-^)10 (a;f-4a;|) 

27Q^tjf I 27q'^iiujf 8l£VM^__|_ 

16/^6 (a;2-4a;|) 4(1-^)1° (tjf-4a;2) 8(l-^i)i" (a;2-4w2) 

27g2^3^3 _ 27g2j^4^3 27gQV(l-M)2x/M^c^f 



4(l-At)V'^('^?-4'^3) 8(l-Ai)V^('^f-4t^3) 16{l-At)i0(cj2- 

9(52ci;2aj3 9q^^uj\uJ3 27q'^ fi'^ui'fujs 

«/ 9 J 9\ J /-. \ 1 n / 9 ^ 9\ r./-. \in/ 9 ^ 



16M«(a;f-4a;§) ^ 4(1-^)1° (aif-4a;f) 8(1-^)1° (a;f-4^ + 

9g2/j3a;2^^3 Qq'^n^ujfujs 9qQy/ (l-fi)^ ^/Jflu)ju3 

4(1-At)l0(w2_4t^2^ 16(1-^)10 (t^?-4w2^ 8(1-m)6m*(wi-4w2) 

9gQ v^cjs'lj3 ggQy/ (1-^^)2 ^^t^f^a , 9g2cji(^g 

4(1-^)^3(0.2-4^32) + 8(l-/.)6M2(^2_4a;32) 2(l-M)io(a;2-4a;|)^ 
9Q^uJiuj^ 18g2^[^]^tj2 27g2/i2(^ji^2 

2m6(^2_4,.32) ~ (1-/.)10(^2_4^2) + (l-;,)10(a,2 4^2y- 

18g^M^a;ia)§ 9g^/j4a;it^g _ 9gQ-^(l-/j)2y^a;ia;2 

(1-A')^ °(^^4^ gj 2(1-m)1° (0^2 4^32) (l-M)VH'^f -40^32) 

18qQy/{l~f,)2y/J^uiiuj^ _ 9gQ^(l-^)2^/^a;icj2 
(1-^)6^.3(0.2-4^^2) (l-^)V2(c.2-4a;2) + 



A 



Slq^uJi I 81Q2a;^ _ 81q^ fiu^i , 81g2^2^^ Sk/Qy/ {l-fi)- y/j^ui 
4(l-/i)io 4/i8 2(l-/i)io 4(l-/i)io 2(l-/i)6^5 



81gQV(l-^)2.^a;i _ 675gv/(l-At)2a,2 675gV(l-M)2^a;f 

2(l-//)6/i4 32(l-/i)8 ^ 32(l-/i)S 

225QV^a>f _ 135g2t^2 ^ 135Q2^2 135g2^a,2 
16//7 16(l-/i)iOcJ3 16^3(^3 + 8(1-At)l0a;3 

135g2;,2^2 ^ 135gQ^(l-/.)2y^a;2 135gQ^(l-M)2A/^a;f 
16(l-/i)i0w3 8(l-/i)6ft''''a;3 + 8{l-fj.)'^ fj.'^a}3 

135g2^(l-^)2^(-l+^)2a,3 _ 135g2^(l-;,)2^(-l+^)2a,3 

4( l-,»)6( - l+M)6 4(1-m)8(-1+m)4 

135g2 ^(1-^)2 y(-l+^)2t^3 ^ 45Q2t^3 135Q2cj3 _90Qfcj3__|_ 

8(1-m)^(-1+/^)* 2(-l+/i)V* 4(-l+^)2^« + (-1+m)V^~^ 

135Q2^3 135Q2t^3 135Q2a;3 90Q2cj3 45Q^W3_ , 

2(-1+m)V^ ( -1+A')V ^ 4(-1+m)^ m'^ + ( _ l+^)4^5 2(-H-m)V''"^ 
135g2^(l-^)2^(_l+^)2^^3 135g2 ^(1-^)2 ^(-i4-^)2^t^3 ^ 
2(1-//)6(-1+;u)6 + 2(1-/^)8 (-1+m)4 + 
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135qV(l-^)V(-l+M)2//u>3 _ 135qV(l-At)V(-l+M)V^u;3 
8 (1-m)7 ( - 1+^)4 4 (l-/x)6( -l+/x)6 

135qV(l-At)V(-l+/x)^M'c^3 45gQV(l-M)^V^a;3 

4(l-^)8(_i+^)4 + 2(1-m)«(-1+m)* + 

135gQi/(l-M)^\/M^c^3 45gQi/(l-M)2y^a,3 

4(1- ^)8 + 2(1-m)6 (-1+|u)V 6 + 

135gQA/(-l+M)^A/M^'^3 _ 1359Qv/(l-At)2y^a)3 

4( -i+m)V ^ ^ 

675gQ^(l-M)^V^u'3 , 135gQ^(l-M)^ ^^^^^3 , 

2(l-M )'^(-l+A t)V* 4(l-/. )«(-l+^Ip 7? ^ 

135gQ^(l-At)2y^c^3 135gQ^(-l+/.)2y;i?a;3 

8(i-^)r(_i+^)2^4 + 4( -l+At)« M'* 

450gQ^(l-At)2^a;3 _ 135qQ^(l-M)^ A/M^a'3 
(l_^)6(_l+^)4^3 (l_^)8(_l+^)2^3 

8(1-M f (-I+m )^'^ 2(l-At )6(-l+M )V2 + 

405qQv'(l-M)2\/A^^3 405gQ^(l-M)2A/^a'3 

2(1-m )»(-1+7 0^ ^ 8(l-At )7(-l+A t)V2 
135qQ.^(l-M)2Y^u)3 _ 135gQ-^(l-/j)2 Y^a;3 

(-i+ /iF7r~ 

135gQ^(l-^t)2./^a)3 405gy^(l-At)2a;ia;3 _ 4:05q^/(l-t^)^ /lUJicJs 

8(1-m)7(-1+m)V 4(1-m)8 4(l-/i)8 + 

135Qv^a>iai3 , 135g^(-l+M)^c^g 45Qv^ai§ _ 225Qv^a)f 

2^7 + 4{-l+At)7 + 2(-1+m)V^ 2(-1+m)V® 

225Qy^cj| 22hQ^/^uil 22bQy/Jfii^l AbQ^/^^l 

T^UWJ? (-l+At)'V^ 2(-1+m)V3 ~ 2PT+^I)V 

75q2a;§ , 25Q2a;| 150Q^aj§ , __375Q\il__ 

(-1+Aj)8^i + (-1+^^)8^-^ 

375gQ^(-l+At)2A/A^t^g 175gQ^(-l+A')^\/M^^| 

2(-l+^t)V'^i 

125gQ.^(-l+/i)2v^c^| 405g2t^; 405g2cjf 

2(-H-M)**/i^wi 8(l-Aj)l2(a;2-4w2^ + 16(1-^4)" (w2-4a;|) 

AObQ'^Lj'f 405g2^i^3 1215q'^fjw'f 

8//8(a;f-4t^|) ~ 2(l-,i)i2 (a;2-4a;|) 16(1-/^)" (a;2-4a;2 j + 

1215g2^2^3 1215g2/f2t^3 405gV^wf_^_ 

4(1-^)12 (a;2 -40,32) + 16(1-^)11 (a;2-4wi) 2(1-^)1 2 (0,2 -4a ,2y 

405g2^3^3 405g2^4^3 ^ 405gQ^(l-^)2 y^Ljf 

16(1-m) "H-^ "^ 8(1-m)^H "i-^'^ 8(l-M)«M«('^f-4c^i) 

405gQ.^(l-/^)2Y^a)3 _ 405gQ-^(l-;u)2y^a.f 

4(l-/.) V5(a;2- 4c.2) 8(1-m) V4(c.!- 4^3) 

405gQ^(l-At)2^/^a;f _ 405gQ^(l-M)2 y/^c^f 

16(l-M)V(^-4a;2) 8(l-/.)VH-4'"3) 

405gQ-y/(l-M)2-v/M2a;f 405gQ^(l-M)2 y^^tjf 
4(l-/.) 8^3(^2_ 4^2j + 16(l-;.)7,,3(<^2_4^2) 

405gQA/(l-/^)^A/A^'^f 135g2i^2t^3 135g2a;2t^3 

8(l-M)V2('^f-4a;32) + 8(l-/.)i2 (c,2 -40.32) + 16(1-^)" (^2-40,32) ^ 
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8M**(t^?-4w2^ ~ 2(l-At)i2(w2-4w2) ~ 16(1-^)11 (a;2-4cj2) + 
405(j2/i2tj2j^3 405l}2;^2^2^g 135^2^3(^2^^ 

4(l-At)l2 (a;2-4a;2) 16(l-/i)ii (cj2^4w.2) ^ 2(1-^1)12 (a;2-4cj2j ~ 

135q2,u3^2^3 135g2/.4^2^3 ^ 135ggi/(l-At)2^c^2^3 

16(l-M)ii(wf-4a)f) 8(l-M)i2(w2-4a;f) 8(l-/i)6/i6(^^2_4^2j + 

135(?Q^(1-m)^a/a^-?'^3 135qQ^(l-M)^V^a^?u;3 



4(1-^)6m"(-?- 


-4-1) 


8(1-m)V*(-?- 


-4-i) 


135gQi/(l-M)2y 


/i^wf(i;3 


135g(3V(l-M)^V 


'/i2t^,2|^3 




-40,2) 


8(l-/.)6/.4(u,2- 




135gQi/(l-^)2y 




135gQ^(l-/.)2y 





135qQ^y (1-^)2 yVu)f c^3 135g2a,ia;2 135g^-3-3_^ 

8(l-M)V2(c^f-4a;i) (1-/.)12 (a;2-4u,2) 2(1-m)11 (a;2-4a;2) ' 

135Q2(j-^(j2 540g2^ci;itJ3 405<j2/ia;itj2 

/.8(a,2_4^2>) + (1_^)12(^2_4^2>) + 2(l-M)ll(c.f-4c^32)~ 

810g2j:^2^^^^ 405g2^2^^^(^2 _540£V-i-3__|_ 

(1-^.)12(,,2_4,^2) 2(l-M)"(c.2-4a;2) + (1-;.)12(<^2_4,^2^ + 



135g2^3^^^2 135g2^4(^,^^2 135(/Qa/ (1-^1)2 y^a;iLj| 

2(1-^)11 (a;2-4c.2) (1-;.)12 (t.2_4<^2>) + (l-^)VB(-?-4a;2) 



(l-/.)6^5(^2_ 


-4-i) 


+ 


(1-M)8/.4(c.2- 


-4-i) 


135gQ^(l-;t)2^ 


/ 9" 2 


+ 


135gQ^(l-,t)2^ 




2(l-,i)7,,4(,^2_ 


-4-1) 




-4^2) 


27(Wy/(l-/,)-\ 


/ //2ljici;3 









(1-/. )V3(,^2 _4^2-) 2(l-/x)Vn-f-4'^3, 
135gQ^ (1-/^)^ V^-1-3 , 270g2cj2a;3 270Q2c^2^3 

{l-,0V2(-?-4c.2) + (l-,.)io(_4c.2+a.2) + p«(-4c.i+c.2)' 



51()(/2/i^2^3 270(i2/i2^2_^3 540(/(^ (!-/')- yV^-i-a 

(l-/.)iu(-4a^cu2) + (1-^)10 (-4a;2+a;2) ^ (1-^) V« (-4o.f +0,2) 

540gQ^ (l-^t)2y^a;2(^3 54g2cjia;2 54Q2a;-^^2 

(l-/i)'V'^(-4-?+-|) ^ (l-,.)l"(-4c.2+a;32) + -'!-?+-') 



108(/2/;_jicj2 51g-/.t2cj^^2 108g(^y/ (l-/;)2 ^/^co'it^g 
(1-^)1" (-4cu2+ cu2) (1-^)10 (-4c^2 +^2) (1-^)6^5 (_4<^2 +^2) 
108gQ^(l~/j)2y^a;m2 54g2a;g 

{l-^,y'^,*{-4ul+u|) (i-A*)"' (-4^^2+^.32) 

54Q2,^3 108g2^t^3 54q2^2^3 

tia{^4ujl+io'i) + (l-/i)l«(-4a;2+a;2) ~ (1-^)10 ( -4<^2_^u)2') " 
108gQi/(l-M)2y^c^3 108gQ^(l-/.)2 
(1-m)VH-4^?+-3) (l-/x)V^(-4a;f+^|) 

when /i -> 0.00025, A 0.00025, g 0.025 andg->0.00025 
D2 = -3.59996 X 10^5 (-30720.2 + ui) coi - + 
(-1.12895 X 10^° + 1.72798 x lO^V) + (^■^^^^^^-^°"-^-^^^^x^°--0- 
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j2(2.21184xl02°tJl+5.5296xl0lSaJ3) 

-i.uj'l+l.uil ' -0.25wf+l.aj| 



(-4.31996 X lOiSwi -4.31996 x lO^^ujs) 



when II 0.00025, A 0.00025, q 0.025, Q->0.00025andw3 ^ 1 Using all the values of 
coefficients and algebraic manipulations we get the following expression for D2'- 
D2 = 5.76096 X 10^^ - M><io^ _ 1.92 x 10^^ui+ 



R v 1 nlO, ,2 I 1.152xl0l5a;i 1.44xl0i^^f 4.32xl0i^a;^ , 

6.b X iu + _4+^^ z:^^ _4+^-2 + 

0.00025 (^-2.30423 x 10^^ + ^-024x10^^ ^ 8.30131 x lO^^^i- 
3.45744 X lO^^u^ - ^-^Y'f.'"'' + ^^^ilM^iofV + 

J- 1— 4cjf 1— 4a;f 



1 ^"'l 

3 



1. 10592x10^-* g;-' _ 5.5296 xlOj'»a;i , 6.912xl0^-^^f 2.0736X lO^^cj^ 

We have computed the value of D2 numerically for various values of parameters when 
/i = 0.0025 and U3 = l.The graphs are plotted D2 versus ui. In figures [H plot I, II, III 
represent the respective values of q2 = 0.25, 0.50, and 0.75 respectively when qi = 0.25 fixed 
and the vertical line shown in each graph is an asymptote. Effect of A2 in which first three 
curves are for A2 = 0.0025 and second three curves belong to A2 = 0.0050. Similarly we 
have also obtained the effect of qi in figure [2] in which plot I and II correspond to qi = 0.50 
and 0.75 respectively when q2 = 0.25 fixed. Here ffist two curves are plotted for A2 = .0025 
and second two curves for A2 = .0050. 




Fig. 1.— D2 when/i ^ 0.00025, A ^ 0.00025, ^ 0.25,0.50, and 0.75, Q ^ 
0.00025 and cug ^ 1 



Dz I Dz n 




Fig. 2.— when 12 0.00025, A -> 0.00025, qi 0.50 and 0.75, Q 0.00025 and 003^1 

3. Conclusion 

It is evident from all the above figures that the curves are in rectangular hyperbolic 
forms with singularity at Ui = .50 and we find that D2 ^ for various values of different 
parameters. Thus we conclude that, according to Arnold's theorm since D2 7^ O,out of plane 
equilibrium point Lq is stable in non-linear sense. 
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